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ABSTRACT. The diversity (Shannon- Weaver) of butterflies throughout the urbanized 
area of Porto Alegre, Brazil, was analyzed using 109 sampling areas within three char- 
acteristic zones of urbanization: buildings (B), houses and buildings (HB), and houses (H). 
Highest diversity was found in the periphery of the houses zone (H) one to two kilometers 
beyond the perimeter of the houses and buildings zone (HB). From zone H to zone HB 
we observed a significant decrease in diversity and a small overlap in community com- 
position (Renkonen’s PS). From zone HB to a central zone of buildings (B), there was a 
relatively small change in community composition demonstrated by statistically similar 
diversity indices and a high species similarity. These findings suggest the existence of two 
macrohabitats for butterflies in the city: 1) B + HB and 2) H. Samples from areas within 
the same urbanized zone showed the highest degree of similarity. Diversity decreased in 
the B + HB macrohabitat mainly owing to a reduction in species richness (S). Biotic and 
abiotic factors that may be involved in this reduction are discussed. For late spring and 
midsummer samples, nearly 50% of the variation in diversity was explained by vegetation 
cover and distance from the city center. This value rises to 63% for the total diversity 
and 70% for the log transformation of species richness. These high proportions emphasize 
the importance of regional urban environmental conditions for butterfly diversity. In the 
winter, only vegetation cover presented a partial regression coefficient that was significant, 
accounting for less than 20% of the variation in diversity. Also, a scattered distribution 
of areas with high butterfly diversity associated with high vegetation cover was observed 
during the winter. 

Additional key words: urban ecology, urban Lepidoptera, insect diversity, man-made 

habitats, community structure. 



The structure and diversity of biotic communities within urban en- 
vironments are important for theoretical and practical reasons. Urban 
biotas can be studied from a genetic and evolutionary perspective, as 
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exemplified by Bishop and Cook (1981), or from an ecological per- 
spective as demonstrated by Ruszczyk (1987). The ecological conse- 
quences of urbanization for particular groups of plants and animals can 
indicate the degree of disturbance of such environments and may be 
useful in developing strategies for conservation. 

In a previous study, Ruszczyk (1987) presented maps of the distri- 
butions of 29 species of butterflies within the urbanized area of Porto 
Alegre, Brazil. Species exhibited variable rates of decline toward the 
highly-developed urban center. The border between a predominantly 
house-occupied zone (H) and a zone of houses and buildings in equal 
proportions (HB) was found to be the main area of transition for the 
urban fauna, acting as an ecological barrier for species typically asso- 
ciated with woods or natural fields. Species that are associated with 
open areas, that are highly vagile, and that have larvae that utilize both 
native and exotic cultivated plants were dominant in the city. Distance 
of the sampling areas from the center of the city was found to be a 
better predictor of butterfly numbers than average elevation or vege- 
tation cover. 

In this study we further analyzed butterfly diversity for 109 sampling 
points. We also investigated factors influencing spatial patterns of di- 
versity and species richness within the city as well as the similarities in 
the structure of butterfly communities in regions within different levels 
of urbanization. 



Materials and Methods 

Study area. Porto Alegre is a large urban area in Rio Grande do Sul 
in southern Brazil (30°02'S 51°14'W; 1,000,000 inhabitants). It has a 
temperate-subtropical climate with high humidity and moderately high 
temperatures in the summer. Mean annual temperature is 13.8°C, and 
average annual rainfall is 1322 mm. Three characteristic zones of ur- 
banization were identified (Ruszczyk 1987): a buildings zone (B) with 
buildings more than four stories high and vegetation cover below 20%; 
a houses-and-buildings zone (HB) with equal proportions of lower build- 
ings and houses and vegetation cover between 20 and 40%; and a houses 
zone (H) with mostly houses but including open areas within the city, 
and vegetation above 40% (Fig. la). The urbanized zones illustrated in 
Fig. 1A were simplified by drawing tangential lines to the borders of 
the different urbanized zones (Fig. lb). 

Data collection. A 1 km grid was superimposed on the map of the 
urbanized zones resulting in 109 contiguous sample areas within the 
city. Sampling areas (SAs) were arbitrarily delineated as 600 m diameter 
circles (Fig. lb). SAs were surveyed for butterflies during three sampling 
periods: November-December 1980, March- April 1981, and June-July 
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Fig. 1. Urbanized zones of Porto Alegre in 1981 (a), and distribution of sampling 
areas in a simplified map (b). The triangle in Fig. lb indicates the center of the buildings 
zone. 1, buildings zone; 2, houses-and-buildings zone; 3, houses zone; 4, marshes; 5, 
extraurban area. 



1981. Each SA was sampled sequentially, five SAs per day between 
1000-1600 h. SAs were censused by walking continually along the streets 
and recording the number of each butterfly species observed during a 
45 minute period (see Ruszczyk 1987 for further details of the sampling 
program and study area). The distance between each SA and the SA 
at the center of the building zone (marked with a triangle in Fig. lb) 
was considered the “distance from the city center/' The average ele- 
vation of each SA was calculated as the arithmetic mean of its highest 
and lowest points. 

Data analysis. The Shannon- Weaver index (Margalef 1958, Lloyd & 
Ghellardi 1964, Pielou 1966) was used to calculate diversity (H 1 ) for 
each SA. The SAs were grouped in three sets related to the three 
urbanized zones. The differences between the calculated indices for 
these sets were compared using a t - test modification proposed by Poole 
(1974) for evaluating diversity calculations. In addition, a one-way 
analysis of variance was applied to the three sampling periods (Nov- 
Dec/80; Mar-Apr/81; Jun-Jul/81) disregarding the zones, to test the 
effects of seasonality on mean diversity. The degree of similarity among 
samples was measured using Renkonen's Percentage of Similarity (PS). 
For this analysis, SAs were combined into 11 regions with similar area 
within a single urbanized zone. This reduced the area matrix from 109 
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Fig. 2. Diversity index (Shannon- Weaver) of butterfly communities in the urbanized 
zones of Porto Alegre. Clear areas indicate less than 12 individuals recorded; bold lines 
separate urbanization zones (see Fig. lb). 



sampling points to 11 contiguous areas and permitted detection of 
patterns among city regions rather than local faunal similarities. Group- 
ing of regions in the principal matrix followed the simple average 
method (Sneath & Sokal 1973). The relative influence of 1) percentage 
of area covered by vegetation, 2) distance from the city center, and 3) 
mean elevation of each SA on butterfly diversity and on the total number 
of species was calculated using multiple regression. The explained vari- 
ation (R 2 ) of the dependent variables was partitioned into components 
attributed to each independent variable following the standard regres- 
sion method (Kim & Kohout 1975). 



Results 

Highest diversity values were found in zone H, typically one to two 
km beyond the perimeter of zone HB (Fig. 2). In late spring and 
midsummer (Figs. 2a & 2b), patterns of diversity were strongly cor- 
related with urbanized zones — zone H typically had diversity values 
greater than H 1 = 2.4 (Fig. 2, diversity class no. 3) and zone HB had 
lower values. During winter (Fig. 2c) there was greater variability in 
the diversity indices. However, the peripheral area of the houses zone 
(H) continued to support greater diversity. 

Papilionini and Heliconiini, two tribes that are abundant in the city, 
showed the same tendencies as described above (Fig. 3). A considerable 
decrease in the number of species in these two groups was observed at 
the border between zones H and HB. In zone B the number of species 
was fairly constant — one or two. The highest number of species was 
recorded in the periphery of zone H. In the winter, the number of 
species decreased greatly, especially in zones B and HB (Fig. 3). 

Multiple regression of diversity (using pooled seasonal data) with the 
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Fig. 3. Number of species of Papilionini (a) and Heliconiini (b) in the urbanized 
zones of Porto Alegre. Each three-column group indicates the number of species recorded 
(from left to right) in late spring 1980, midsummer 1981 and winter 1981. 



number of species and species evenness showed a greater interaction 
with the number of species in the explained variance of diversity (R 2 
= 0.90) than with the species evenness (R 2 = 0.10). The standardized 
regression coefficient (variables transformed to have unit variance al- 
lowing the comparison of variables measured in different units) of the 
number of species (0.98) was nearly three times greater than for even- 
ness (0.35), indicating that diversity decreased in zones B and HB mainly 
due to the reduction in the number of species in these zones. Diversity 
indices for the entire zone H (summation of SA data) were significantly 
higher than for zones HB and B. On the other hand, the indices for 
zone HB and B were not statistically different (Table 1). 

Table 1. Diversity of butterflies in three urbanized zones of Porto Alegre. The lines 
indicate indices that were statistically similar (P < 0.05). Differences between indices 
were compared using a t-test modification proposed by Poole (1974). 







Urbanized zone 




Sample 


Buildings 


Houses-and- 

buildings 


Houses 


Late spring, 1980 


2.5414 


3.0304 


3.4785 


Midsummer, 1981 


2.6344 


3.0145 


3.3108 


Total 


2.9277 


3.2251 


3.4873 
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Fig. 4. Similarity (Renkonen’s index) in the structure of butterfly communities of 
regions within the city of Porto Alegre. I to VII, houses zone; VIII to X, houses-and- 
buildings zone; XI, buildings zone. 



Figure 4 illustrates the results of PS calculations. Two distinct clusters 
of regions are evident: 1) regions I to VIII, all but the last of which are 
located in zone H, with a high internal PS; and 2) regions IX and X in 
zone HB and region XI in zone B, with a comparatively low internal PS. 

Standardized regression coefficients of distance from the city center 
and vegetation cover with diversity and the log number of species, were 
similar; but distance had a slightly greater value (Table 2). In midsum- 
mer, distance was almost twice that of vegetation cover. In the winter, 
vegetation cover showed a significant partial regression coefficient with 
diversity, though accounting for less than 20% of the explained variation 
in diversity. In late spring and midsummer, about 50% of the variation 
in diversity was explained by vegetation cover and distance from the 
city center. This value increased to 63% for the total diversity and 70% 
for the log number of species. These high proportions demonstrate that 
such variables are important determinants of butterfly community 
structure in Porto Alegre. However, the contribution of each variable 
alone was, in general, less than 11%. More than 30% of the variation 
was due to the interaction of vegetation cover and distance, whose 
effects on butterfly diversity and species number was neither indepen- 
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Table 2. Decomposition of the explained variation in butterfly diversity (Shannon- 
Weaver index) in three sampling periods, and total number of species (In S) into com- 
ponents attributed to the independent variables percentage of area covered by vegetation 
of the sampled area (Xi) and distance from the city center (X 2 ). Standardized regression 
coefficients in parentheses. Angular transformation was used for Xi and logarithmic for X 2 . 



Dependent 

variable 

(Y) 


Proportion 
of variation 
explained by 
Xi and X 2 
(R2) 


Increment due 
to vegetation 
cover 

(Xi) 


Increment due to 
distance from the 
city center 

(X 2 ) 


Not attributed 
to either Xi 
or X 2 alone 


Diversity 










Late Spring 


0.456 


0.058 


0.075 


0.323 


1980 




(0.341) 


(0.390) 




Midsummer 


0.489 


0.043 


0.108 


0.338 


1981 




(0.289) 


(0.464) 




Winter 


0.171 


0.171 


— 


— 


1981 




(0.413) 






Total 


0.630 


0.078 


0.107 


0.445 


1980-1981 




(0.394) 


(0.465) 




Log no species 


0.703 


0.096 


0.109 


0.498 






(0.440) 


(0.467) 





( — ), partial regression coefficient not significant after the inclusion of the variable in the multiple regression equation. 



dent nor additive (Table 2). The biological mechanism underlying the 
strong association between these two variables affecting butterfly di- 
versity and species number cannot be inferred from these data alone. 
Average elevation was excluded from this analysis because in the pres- 
ence of the other variables its partial regression coefficient was not 
significant. 

A marked seasonality in diversity was observed in Porto Alegre (Fig. 
5), with mean diversity of each sampling period significantly different 
even when the zone data were pooled (Table 3). 

Discussion 

The buildings zone (B) did not form a subset isolated from the houses- 
and-buildings zone (HB) (Fig. 4), and the diversity indices for these 
zones were statistically similar (Table 1) revealing their similar butterfly 

WINTER 1981 LATE SPRING1980 MIDSUMMER 1981 



H=2.242 I-L2.354 H=2.532 




-1 ! 1 ! — f— 1 1 — 1 1 — I H 

2.1 2.2 2.3 2.4 2.5 2.6 H 1 

Fig. 5. Confidence intervals (P = 0.95) of the means of butterfly diversity in urban 
Porto Alegre during three sampling periods. Number of sampled areas between paren- 
theses. 
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Table 3. Results of the one-way ANOVA for mean diversity in three sampling periods, 
independently of urbanized zone. Hi = Late Spring 1980; H 2 — Midsummer 1981, H 3 
= Winter 1981. 



Source 


ss 


DF 


MS 


F 


P 


Between 


3.7050 


2 


1.8525 


15.47 


< 0.001 


(H r H 2 ) vs. H 3 


2.1274 


1 


2.1274 


17.77 


< 0.001 


Hi vs. H 2 


1.5776 


1 


1.5776 


13.18 


< 0.001 


Within 


31.9676 


267 


0.1197 






Total 


35.673 


269 









community structure. Zone B can be seen as a portion of the zone B 
+ HB macrohabitat where environmental conditions are more harsh 
for butterflies, but without divergence of the typical community of 
butterflies of zone HB. In this community (B), several species of but- 
terflies that are relatively abundant in zone H, are scarce or absent. 
These include woodland species, such as Heliconius erato phyllis (Fa- 
bricius), Eunica margarita (Godart), Adelpha spp. (Nymphalidae), Bat- 
tus spp., Parides spp. (Papilionidae); Nymphalidae that feed on fruit 
and sap ( Hamadryas spp., Anaea spp.); species characteristic of fields, 
such as Colias lesbia pyrrhothea (Huebner), Eurema spp. (Pieridae), 
Euryades corethrus (Boisduval) (Papilionidae), Junonia evarete (Cra- 
mer), Vanessa spp. (Nymphalidae); and some eurytopic species, such 
as Dryas iulia (Fabricius), Anartia amathea (Eschscholtz) (Nymphal- 
idae), Papilio hectorides Esper, and Papilio astyalus Latreille (Papili- 
onidae). All species found in zones B and HB also were observed in 
zone H. These species [e.g., Papilio scamander Boisduval, Papilo an- 
chisiades capys Huebner (Papilionidae), Ascia monuste orseis (La- 
treille), Tatochila autodice (Huebner) , Phoebis philea (Johansson) (Pier- 
idae), and Dryas iulia] are the most widespread in the city, attaining 
high densities in all urbanized zones (Ruszczyk 1987). The impoverish- 
ment of the butterfly community in the B + HB urbanized zone is 
likely the result of the considerable environmental disturbance in this 
area as compared to zone H. Abiotic and biotic factors are harsher in 
zone B + HB than in zone H. Abiotic factors include the following: a) 
dry and strongly illuminated habitat (e.g., all watercourses are chan- 
nelized, there are few shaded areas, and there is high sunlight pene- 
tration to the ground); b) greater air pollution owing to traffic; c) habitat 
disturbance from intense human movement and traffic; and d) streets 
and sidewalks completely paved and a large percentage of the area 
occupied by buildings, diminishing the resources at the soil surface. 
Biotic factors include the following: a) lower percentage of area covered 
by vegetation which acts to decrease diversity within these zones by 
lowering primary productivity (Connell & Orias 1964), and the accen- 
tuated fragmentation of the vegetation probably reduces colonization 
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and dispersion of butterflies; b) relative homogeneity of vegetation (in 
zones B and HB relicts of native vegetation were not observed, and 
many plants that are common in outlying areas were very scarce, giving 
a qualitative decrease in nectar sources and potential food plants for 
larvae); and c) smaller contribution of elements of the extraurban fauna 
(zones B and HB probably have small participation of the transitory 
species from peripheral areas than zone H, which is in direct contact 
with remnants of natural habitats). 

The formation of butterfly diversity gradients in the city is in contrast 
to the results obtained in studies of soil arthropods, which seem to 
respond more to local (soil) variables than to urban environmental 
gradients (Kuhnelt 1955, Topp 1972, Lussenhop 1973, Maurer 1974). 
More recently, Klausnitzer and Richter (1983) demonstrated the pres- 
ence of an urban gradient for carabids in the city of Leipzig, Germany. 
As for butterflies, the distance from the center of the city showed a 
greater influence on diversity than vegetation cover or mean elevation, 
emphasizing the importance of macrohabitat conditions for these in- 
sects. 

Distance correlated well with butterfly abundance (Ruszczyk 1987) 
and diversity (present paper) in Porto Alegre, Brazil, probably because 
many parameters that are important for butterflies are radially dis- 
persed in the city due to the radial pattern of the urbanized zones. 

The non-significance of the partial regression coefficient of distance 
from the city center during winter, and the island character of the class 
distribution of diversity (Fig. 2c), suggest that winter butterfly diversity 
depends on qualitative variables such as presence of habitat refugia. 
When compared to other zones, zone H showed greater possibilities for 
the presence of such refugium. In the winter, a higher diversity was 
scattered among 26 SAs, all but one of which was situated in zone H; 
19 of these SAs possess vegetation cover greater than 45%. This provides 
evidence of the biotic value of urban vegetation, and suggests that 
fragments of habitat within the urbanized areas, especially urban forest 
fragments (Rodrigues et al. 1992), may perform a vital role in main- 
taining local biodiversity. 
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